Anomalous Binder Cumulant and Lack of Self-averageness 
in Systems with Quenched Disorder 
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The Binder cumulant (BC) has been widely used for locating the phase transition point accurately 
in systems with thermal noise. In systems with quenched disorder, the BC may show subtle finite- 
size effects due to large sample-to-sample fluctuations. We study the globally coupled Kuramoto 
model of interacting limit-cycle oscillators with random natural frequencies and find an anomalous 
dip in the BC near the transition. We show that the dip is related to non-self-averageness of the order 
parameter at the transition. Alternative definitions of the BC, which do not show any anomalous 
behavior regardless of the existence of non-self-averageness, are proposed. 
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The characterization of phase transitions relies mainly 
on the singularity structure of physical quantities at the 
transition, which can be quantified by critical exponent 
values. In numerical efforts, the accuracy of the esti- 
mated exponents heavily depends on the precision of lo- 
cating the phase transition point. In the case of most 
thermal systems, the Binder cumulant (BC) is widely 
believed to provide one of the most accurate tools for es- 
timating the transition point [H, 0, 0|. The critical BC 
value at the transition is also believed to be universal, 
even though there is still controversy over its universal- 
ity 1. 

In some complex systems [5| , the BC shows an anoma- 
lous negative dip in finite systems, which represents a 
rugged landscape (multi-peak structure) in the probabil- 
ity distribution function (PDF) of the order parameter. 
Great care is required in analyzing numerical data to see 
whether the dip will vanish in the thermodynamic limit. 
If it does, the negative dips in the finite systems can be 
attributed to long-living metastable states. Otherwise, a 
nonvanishing negative dip usually implies that the tran- 
sition is not continuous, but is of the first order. 

In systems with quenched disorder, the disorder fluc- 
tuation may also generate an anomalous negative dip in 
the conventional BC, which is defined as the ratio of the 
disorder-averaged moments of the order parameter. In 
this case, the negative dip may be related to the non-self- 
averageness (NSA) of the order parameter, which usually 
implies an extended and/ or multi-peak structure in the 
disorder-averaged PDF [6[ . 

We consider a typical nonequilibrium dynamical sys- 
tem with quenched disorder, such as the Kuramoto model 
of interacting limit-cycle oscillators with random natural 



frequencies jjj]. The dynamic synchronization transition 
is dominated by space-time fluctuations of the order pa- 
rameter. The quenched disorder is, by definition, per- 
fectly correlated in the time direction, so it may generate 
strong disorder fluctuations similar to quantum systems 
with random defects [8(. In fact, we recently showed that 
the disorder fluctuation was anomalously strong near the 
synchronization transition [^, [13] • 

We take the globally coupled Kuramoto model, which 
can be solved analytically to some extent. The model is 
defined by the set of equations of motion 



dt 



i 



(1) 



where ipt represents the phase of the ith limit-cycle os- 
cillator (i = 1, 2, • ■ • , N). The first term u>i on the right- 
hand side denotes the natural frequency of the ith oscil- 
lator, where u>i is assumed to be randomly distributed ac- 
cording to the Gaussian distribution function g(u>) char- 
acterized by the correlation (LUiUj) — 2o~Sij and zero 
mean ((wi) — 0). 

We note that the natural frequency u>i plays the role 
of "quenched disorder". The second term of Eq. JT]) rep- 
resents global (all-to-all) coupling with equal coupling 
strength K/N . The sine coupling form is the most gen- 
eral representation of the coupling in the lowest order of 
the complex Ginzburg-Landau (CGL) description and 
its periodic nature is generic in limit-cycle oscillator sys- 
tems. We consider the ferromagnetic coupling (K > 0), 
so the neighboring oscillators favor their phase difference 
being minimized. The scattered natural frequencies and 
the coupling of the oscillators compete with each other. 
When the coupling becomes strong enough to overcome 



2 



the dispersion of natural frequencies, macroscopic regions 
in which the oscillators are synchronized by sharing a 
coupling-modified common frequency Q — may emerge. 

Collective phase synchronization is conveniently de- 
scribed by the complex order parameter defined by 
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where the amplitude A measures the phase synchroniza- 
tion and 9 indicates the average phase. When the cou- 
pling is weak, each oscillator tends to evolve with its own 
natural frequency, resulting in the fully random desyn- 
chronized phase (A = 0). As the coupling increases, some 
oscillators with uji w become synchronized, and their 
phases 0, start to show some ordering (A > 0). 

Equation {I} can be simplified to A decoupled equa- 
tions 



d(fj 
dt 



KAmn(ifi - 9), 



(3) 



where A and 9 are to be determined by imposing self- 
consistency. In the steady state (t — > oo), the self- 
consistency equation reads 



A = a(KA) - b{KA) 3 + <D(KAf 



(4) 



with a = (n/2)g(0) and c = -(n/16)g"(0) 7]. This 
equation has a nontrivial solution only when K > K c — 
I/a: 



A ~ (K - K c f 



(5) 



with f3 = 1/2. We note that the exponent (3=1/2 
corresponds to the mean field (MF) value for systems of 
locally coupled oscillators 0. 

Now, we perform numerical integrations of Eq. ([1]) 
by using Heun's method [ll[ for various system sizes of 
A = 200 to 12800. For a given distribution of disorder 
{u>i}, we average over time in the steady state after some 
transient time. After the time average, we also average 
over disorder. Typically, we take the time step St = 0.05, 
the maximum number of time steps N t = 4 x I0 4 , and the 
number of samples N s = 100 ~ 1000. For convenience, 
we set 2cr = 1 (unit variance); then, the corresponding 
critical parameter value is K c = \/&pK = 1.595769.... 

Figure [T] shows the behavior of the phase synchroniza- 
tion order parameter A against the coupling strength K 
for various system sizes A. In the weak coupling region 
(K < 1.6), we find that the order parameter approaches 
zero as A - A" 1 / 2 , which is a characteristic of the fully 
random phase. In the strong coupling region (K > 1.6), 
A saturates to a finite value, indicating a phase transi- 
tion at K w 1.6 in the thermodynamic limit (A — > oo), 
which is consistent with the analytic result. 

To pin down the transition point K c precisely, we use 
the Binder cumulant method [l], 0). The fourth-order 
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FIG. 1: Phase synchronization order parameter A versus the 
coupling strength K for various system sizes N. 



cumulant of the order parameter, the Binder cumulant 
(BC), is defined in thermal systems as 



B • 



1 



(A 4 ) 
3(A 2 ) i 



(0) 



where (• ■ • ) represents the thermal (time) average. In 
systems with quenched disorder, on the other hand, we 
should consider the disorder average besides the ther- 
mal one. We may first consider the BC as the disorder- 
averaged moment ratio [H, 0, EH 



B 



(i) 

A 



[(A 4 )] 
3[<A 2 )F 



(7) 



where [• ■ • ] denotes the disorder average, i.e., the average 
over different realizations of {ui}. 

Figure [2] displays Bj± as a function of the coupling 
strength K for various system sizes N . In the region of 
weak coupling (K — ► 0), we expect the random nature of 
the oscillator phases {<pi} to yield an asymmetric Poisson- 
like probability distribution function (PDF) character- 
ized by P(A) ~ Aexp(— cA 2 ) with a constant c, which 
leads to B^ = 1/3. On the other hand, in the strong- 
coupling region, the PDF becomes a <5-like function with 
a very narrow variance, which leads to B^ =2/3. The 
numerical data in Fig. 2 are consistent with our predic- 
tions. 

However, near the transition, the B 1 ^ shows a big 
anomalous "dip" on the desynchronized side. As the sys- 
tem size increases, the dip develops initially with a broad 
width and then becomes sharper and also deeper. The 
dip's position moves toward the transition point. The 
crossing points seem to nicely converge to the critical 
point K c = However, as the system size increases, 

the presence of the dip starts to hinder us in locating the 
critical point accurately. 

In this Letter, we explain why the dip develops in this 
system and propose alternative definitions of the Binder 
cumulant that do not show any dip in the same system. 
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FIG. 2: Binder cumulant b£ } , denned by Eq. ©, of the 
phase synchronization order parameter A. 
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FIG. 3: Disorder fluctuation A A 2 defined by Eq. (JBJ. 



We measure the disorder (sample-to-sample) fluctuations 
denned as 



A = 



nor 
now 



i. 



(8) 



where O is any observable, such as A and A 2 , in a sys- 
tem. This quantity is positive definite and is supposed 
to vanish in the thermodynamic limit in self-averaging 
systems and to remain finite in non-self-averaging sys- 
tems 0, As one can see in Fig. [3j the disorder 
fluctuation A^2 is quite sizable in the range of K where 
the dip appears (A& shows a similar behavior). In other 
words, the shows a dip where the system is not well 
self-averaged. A careful finite-size analysis on A&2 re- 
veals that it vanishes as ~ TV -1 away from criticality, 
but saturates to a finite value at criticality. The non- 
self- aver ageness at criticality is not surprising because 
the quenched randomness in natural frequencies should 
be relevant at this transition. 

Strong disorder fluctuations may cause non-negligible 
spreading of the effective coupling constants over differ- 
ent realizations of disorder [l3j . Figure 0] shows for 20 
independent samples, the PDF of A just below the tran- 
sition and obtained from the time series of A after the 
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FIG. 4: Probability distribution function (PDF) at K = 1.5 
and N = 12800 in the steady state. Each curve corresponds 
to one of 20 independent samples. 
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FIG. 5: Binder cumulant B^' defined by Eq. ©. Note that 
the dip shown in B'x disappears. 



system had reached the steady state. Indeed, a large part 
of the sample-to-sample variations can be interpreted as 
a shift in the K c of individual samples. The two quan- 
tities [(A 2 )] and [(A 4 )] in Eq. §7§ can be considered 
as the second and the fourth moments of the disorder- 
averaged PDF, which is much broader than the individual 
PDFs near the transition. One can easily see that broad- 
ening yields a larger value for the ratio [(A 4 )] /[(A 2 )] 2 
and, hence, a smaller BC. The effect is particularly pro- 
nounced on the small K side of the transition, where A 
itself is small, in which case a shift in K c has a stronger 
influence on the moments. 

An alternative definition for the Binder cumulant 
for systems with quenched disorder (especially non- 
diminishing disorder fluctuations) isfsLflij] 



D 



(2) 



(A 4 ) 
3(A 2 ) S 



(9) 



We note that the disorder average is performed over the 
ratio of the time-averaged moments. The moment ratio is 
calculated for each sample first and, is then averaged over 
disorder. It is clear that this definition of the Binder cu- 
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FIG. 6: Binder cumulant B^ 3) defined by Eq. flQjl 



mulant should eliminate the most dominant contribution 
from the disorder fluctuations, i.e., the anomaly caused 
by the spreading of the effective coupling constants. This 
definition has been adopted mostly in quantum disorder 
systems, where strong disorder fluctuations are antici- 
pated Q. Figure E] displays B^ A versus K. We note that 
the dip shown in Fig. [2] disappears and that the crossing 

(2) 

points nicely converge to K c , implying that B A should 
serve better for locating the transition point than the 
conventional one, which is confirmed numerically (not 
shown here). 

Yet another definition of the Binder cumulant is 



B 



(3) 



A 



= I - 



[(A 4 )] 
3[(A 2 ) 2 ] 



(10) 



We expect that B\' may also behave smoothly near the 
transition because it does not involve disorder fluctuation 
terms such as [• ■ -] 2 included in B A . Figure [6] displays 
B\ versus K. As expected, we find no anomalous be- 



havior in B A \ We can directly relate B A ' and B 



(i) 



;(3) 



through the disorder fluctuation A&2. 
leads to 



Simple algebra 



B 



(i) 



B 



(3) 



(l-B^)A 



A 2 - 



(11) 



(1) 



As the disorder fluctuation A&2 becomes larger, B A 
shows a bigger dip. This explains quantitatively the 
size and the location of the dip in B A '. The critical 

value of B A (« 2/3) provides additional information 
on the temporal variations of A. One can show that 
3B { 1 ] = 2-[<<5A 2 )]/[(A 2 ) 2 ], where (<5A 2 ) = (A 4 )-(A 2 ) 2 . 
Our numerical result indicates that the relative temporal 
fluctuations are almost negligible even at criticality. In 

(3) 

this case, B A is not practically useful in locating the 
transition point accurately. 

In summary, we studied Binder cumulants in the 
quenched disorder system. For the Kuramoto model, we 
found that the conventionally defined BC shows a big 
anomalous dip near the transition. This dip is shown to 



be directly related to the disorder fluctuation (non-self- 
averageness). Alternative definitions of the BC, which 
did not show any anomalous behavior were proposed and 
may be useful in locating the transition point accurately 
in general systems with quenched disorder. 
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